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Abstract 

We propose a new method for extracting neutron structure functions from inclusive structure 
functions of nuclei, which employs an iterative procedure of solving integral convolution equations. 
Unlike earlier approaches, the new method is applicable to both spin-averaged and spin-dependent 
structure functions. We test the reliability of the method on unpolarized F2 and polarized gi 
structure functions of the deuteron in both the nucleon resonance and deep inelastic regions. The 
new method is able to reproduce known input functions of almost arbitrary shape to very good 
accuracy with only several iterations. 
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I. INTRODUCTION 



Understanding the detailed structure of the nucleon remains one of the central problems 
of the strong nuclear interactions. This is particularly challenging in the so-called transition 
region at momentum scales ~ 1 GeV, where neither perturbative Quantum Chromo dynamics 
(QCD) nor effective hadronic theories provide adequate descriptions of physical observables. 

Over the past few years one of the fascinating developments in the study of this transi- 
tion has been the phenomenon of quark-hadron duality in inclusive electron-nucleon (and 
nucleus) scattering. Here the structure functions in the region dominated by low-lying reso- 
nant excitations of the nucleon are found to closely resemble, on average, the deep inelastic 
structure functions describing the high-energy cross section Since QCD at high energy 
and momentum transfers can be treated perturbatively, but is highly nonperturbative at low 
energies where hadronic degrees of freedom are prominent, this duality provides an intimate 
link between the two regimes. 

Recent experiments have sought to quantify quark-hadron duality by determining its 
flavor, spin, and nuclear dependence, while theoretical endeavors have attempted to un- 
derstand its dynamical origin from a more fundamental basis (for a review see Ref. j^). 
In particular, even from simple quark model arguments one expects intriguingly different 
behaviors of duality for the proton and for the neutron 

While duality for the proton has been tested to rather good accuracy in recent mea- 
surements, for both unpolarized and polarized scattering 0, Isl, there is almost a complete 
absence of analogous empirical information on the neutron. This lack of knowledge has 
prevented the various theoretical models from being adequately tested, and has impeded 
progress in unraveling the microscopic origin of the duality phenomenon. The difficulty 
with obtaining data on neutron structure functions is of course the absence of free neutron 
targets. As a result one often makes use of light nuclei such as deuterium 3 or 3He Q as 
effective neutron targets, assuming that the nuclear corrections are negligible. 

Even when nuclear effects are considered, there exist practical difficulties with extracting 
information on the neutron from nuclear data. Some attempts have been made to obtain 
the spin-averaged structure function from proton and deuterium data in the deep in- 
elastic scattering region, where the exchanged four-momentum transfer squared is large 
(~ 5 GeV^ or greater). A common approach has made use of the so-called smearing factor 
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method |8| , where after an initial guess for one iterates the solution in order to eliminate 
the dependence of the extracted neutron structure function on the starting point. 

In practice the smearing factor method has only been applied to the unpolarized F2 
structure function and only in deep inelastic kinematics 9|. The robustness of this procedure 
is guaranteed only for functions which do not change sign, and for spin-dependent structure 
functions, which can have several zeros, the usual prescription is inadequate. Furthermore, 
in the nucleon resonance region, where there exists non-trivial resonant structure, it is not a 
priori clear whether it is even possible to extract resonance structure that has been smeared 
out by nucleon Fermi motion. 

In this paper we propose a new method in which the nuclear effects are parameterized 
via an additive correction to the free nucleon structure functions. In contrast to the more 
common multiplicative method, which is problematic for structure functions with zeros, 
the new method can be used for functions of almost arbitrary shape, which allows access to 
neutron structure in both the deep inelastic and resonance regions. By iterating the solution, 
the dependence on the initial guess for the neutron structure function is eliminated, and in 
practice a reliable extraction can be achieved after only several iterations. 

In Sec. |TT]we present the formalism for computing nuclear structure functions at finite 
within the nuclear impulse approximation. While the formalism is general and can be applied 
to any nucleus, to illustrate the features of the new extraction method we focus on the specific 
case of the deuteron. In Sec. Illll we present the details of the new method, and discuss other 
methods which have been used to extract neutron structure functions from nuclear data, 
including the smearing factor method, and a direct method of inverting integral equations 



which has previously been used in Refs. [lO|, Our results are presented in Sec. [IV] 



for spin-averaged and spin-dependent structure functions, in both the resonance and deep 
inelastic regions. Using known input functions constructed from resonance and leading twist 
structure function parameterizations, we demonstrate the accuracy of the extraction method 
and provide a detailed discussion of its convergence to the exact results as a function of the 
number of iterations and the first guess in the iteration. Finally, in Sec. |V] we summarize 
our results and preview future applications of the new method. 
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II. NUCLEAR STRUCTURE FUNCTIONS 



The usual framework for computing structure functions of nuclei at large x is the rela- 
tivistic nuclear impulse approximation, in which the lepton probe scatters from the nucleus 
incoherently via the scattering from its bound proton and neutron constituents. In this 
approximation the nuclear structure functions can be written as convolutions of the bound 
nucleon structure functions and nucleon light-cone momentum distributions in the nucleus 
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In particular, for the spin-averaged F2 structure function of a nucleus A we have: 

F,\x, Q') = (/o^/^ ® F|) (x, Q') + [f-/^ ® F^) (x, Q') , (1) 

where x = Q"^ /2Mav is the Bjorken scaling variable (per nucleon). Ma is the nuclear mass 
and V is the energy transfer, and the symbol ® denotes the convolution 



T ® F,-) Q') - / dy f^'\y^ 7) Q') , (2) 



with M the nucleon mass. The function f^^"^ is the light-cone momentum distribution 
of nucleons in the nucleus, and is a function of the light-cone momentum fraction y of 
the nucleus carried by protons (A^ = p) or neutrons (A^ = n), and of the virtual photon 
"velocity" 7 in the target rest frame, 7 = \q\/qQ = (1 + 4M^x^/(5^)^/^. For moderate 



values ~ 1 — 10 GeV^, 7 ranges between unity and ~ 2. As discussed in Ref. [l8|, Il9| . 
taking the full dependence of the smearing function into account is vital for discussing 
nuclear structure functions at large-x or resonance kinematics. 

For the spin-dependent nuclear and (7^ structure functions, one has a set of coupled 
equations involving both the nucleon and g2 structure functions [l^, 

xgfix. Q') = [ft ® ^9^) (^^ Q') + {fT ® ^sf) (x, Q') , 2, J = 1, 2 (3) 

where fl^^^ are the spin-dependent nucleon light-cone momentum distribution functions in 
the nucleus, and a sum over repeated indices j is implied. In contrast to F.^, which receives 
contributions only from the nucleon structure function, the spin-dependent structure 
functions at finite involve also non-diagonal contributions f^^^ and f^/^. (Note that 
both the transverse F^ and longitudinal F^ structure functions individually receive non- 



diagonal contributions, whereas F^ does not 



16j.) In the Bjorken limit, the distribution 



fi2^ vanishes, and the expression for the gi structure function becomes diagonaL Equa- 
tions ([1]) and ([3]) can be viewed as equations in the single independent variable x for fixed 
values of Q^. In the following, for ease of notation we suppress the dependence of the 
structure functions on Q^. 

The light-cone momentum distribution functions in Eqs. ([I]) and ([3]) (also referred to 
as smearing functions) can in general be calculated from nuclear spectral functions which 
account for the ground state wave function of the nucleus and the excitation spectrum of 
the spectator nuclear system, including the continuum spectrum. Since the characteristic 
energies and momenta of the bound nucleons are small compared with the nucleon mass M, 
the unpolarized distribution /o and the spin-dependent diagonal distributions /n and /22 
are sharply peaked about y = 1. 

In this analysis we will focus on the case of the deuteron, for which the smearing functions 



have recently been evaluated in the weak binding approximation [16|, Il8| , including the finite- 
corrections encoded through the dependence on 7. Note that in the isospin symmetric 
limit the proton and neutron distributions in the deuteron are identical, /^^'^ = /"^'^ = /, 
and we shall in the following omit the superscripts on these distributions. 

The unpolarized /o distribution function is given in terms of the deuteron wave function 



Mp) by (i6|,1i7| 



where e = — / (2M), with Sd = —2.2 MeV the deuteron binding energy. The analogous 



spin-dependent light-cone distributions fij are given explicitly in Ref . [18| . For 7 = 1 the /o 
distribution is normalized to the number of protons or neutrons in the deuteron, while /n 
is normalized to the nucleon polarization in the deuteron, 

dyfoivA) = 1 , (5) 



Ma/m 2 

dy fii{y, I) = 1 - - ujd , (6) 



where Ud is the deuteron D-state probability. For the Paris deuteron wave function [20 1 
used here Ud = 5.8%. At finite Q^, or 7 > 1, these normalization conditions are no longer 
satisfied, and the distributions do not have probabilistic interpretations. 





FIG. 1: (Color online) Nucleon distribution functions in the deuteron for 7 = 1 (solid), 1.5 (dotted) 
and 2 (dashed): (a) fo{y,^) distribution for the F2 structure function, (b) /ii(y, 7) distribution 
for the xgf structure function. 



In Fig. [T] we show the /o(y, 7) distribution for the unpolarized structure function 
[igI and the diagonal /ii(y, 7) distribution for the polarized xgf structure function 18|, for 
7 = 1, 1.5 and 2, using the deuteron wave function obtained from the Paris nucleon-nucleon 
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potential [20|. For 7 = 1 the (Q^-independent) distributions are peaked sharply around 
y = 1, and decrease rapidly with increasing |y — 1|, so that by |y — 1| > 0.4 they become 
almost negligible. For larger 7 the distributions become broader, with a larger width and 
smaller height at the peak. For 7 = 2 the height of the peak is approximately half of that 
for 7 = 1. 

The shapes of the unpolarized and polarized distributions are very similar, with ^ 0.5% 
differences between /o and /n at the peak for 7 = 1, and < 2.5% for 7 = 2. While the 
spin-averaged function /o is constrained to be positive, the spin-dependent /n function, 
which involves a difference of distributions of nucleons with spins aligned and anti-aligned 
with that of the deuteron, need not be positive. For values oi y < 0.8 the 7 = 1) 

distribution in fact becomes slightly negative, as is (barely) visible in Fig. [T], although the 
smearing functions here are close to zero. 

Before proceeding to the discussion of extraction methods using these distributions, we 
should note that while the impulse approximation (scattering from individual nucleons in 
the nucleus) provides the main contribution to nuclear deep inelastic scattering, in realis- 
tic calculations of nuclear structure functions other effects are also known to play a role. 
These include nuclear shadowing and meson exchange currents at small x, final state inter- 
actions of the produced hadronic state, relativistic effects, and off-shell corrections to the 
bound nucleon structure functions [l^. Some of these may be formulated within general- 



ized convo 
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utions, either as two-dimensional convolutions with off-shell structure functions 



151], or in terms of exchanged- meson smearing functions 21|, and the techniques 



discussed here may be applicable. Others, such as relativistic corrections, go beyond the 
convolution approximation [l^, and must be included as additive corrections to the con- 
volution. Explicit calculations of final state interactions in the quasi-elastic region have 
suggested that rescattering effects decrease with increasing 22], and in addition partly 
cancel in inclusive inelastic cross sections when summed over several exclusive channels |23| . 

In the present analysis we do not attempt to provide a complete description of nuclear 
structure functions; instead we wish to study the usefulness of the new method of unsmear- 
ing nucleon structure functions within the conventional convolution framework. Once we 
establish the methodology of the new method, additional effects beyond the convolution 
approximation can be considered in actual data analyses. 
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III. EXTRACTION METHODS 



Having outlined the formalism for computing structure functions of nuclei in terms of 
those of nucleons, in this section we review several methods for extracting neutron structure 
functions from nuclear (in practice, deuterium) data, including the new "additive" method 
proposed in this paper. To extract the neutron structure function from proton and nuclear 
data at a given Q^, one first convolutes (or smears) the proton structure function JF^, where 

= F2 or xgi^2, with the appropriate smearing function, 

T^x) = (/ ® j^n (^) , (7) 

where f = fo for the unpolarized F2 structure function, and / = fij for the polarized gi^2 
structure functions. Subtracting the smeared proton J-'^ from the nuclear structure function, 
one obtains an effective smeared neutron structure function 

^"(x) = J^''{x) - ?%x) , (8) 

and then solves the equation 

^"(x) = (/®^")(x) (9) 

for jF"'(x). Note that for a fixed Q^, 7 is a function of x alone, so in practice the smearing 
functions acquire an x dependence. 



A. Direct solution 



Equation ([9]) is a system of so-called Volterra integral equations of the first kind, which 
take the general form 

ry max 

g{x)= / dyKix,y)ziy) , (10) 



where g{x) and K{x, y) (the kernel) are known functions and z is unknown. The general 



theory of Volterra equations is quite extensive, see for example Ref. 2J]. Most Volterra 
equations have no closed-form solution, but numerical solutions for first-kind equations are 
quite simple. Dividing the interval < y < ?/max into a grid of width h hj ya = ah, with 
a = 0,1, . . . , N, and using a quadrature method such as the trapezoidal rule or Simpson's 
rule, one can approximate the integral in Eq. (fTOj) by a discrete sum 

N 

ga = ^I<^ab Zb , (11) 

b=a 



reducing the numerical solution to a problem of matrix inversion: z = K~^g. In fact, 
because of the variable lower limit of integration y = x, the matrix K is upper-triangular, 
and the inversion is almost trivial. This method has been utilized in Refs. 10|, lUl] in a 
similar application. The method fails, though, if K is singular. 

Letting t = x/y and v = x/ymax, Eq- Q can be expressed in the form of Eq. (ITUi) : 

^"(ymaxtO =fdtf (^,7) ^ -^"(t) , (12) 

in which case the kernel is a sum of terms K{v,t) = {ymiaV/t^) /(^max'^/^, 7)- The diagonal 
K{v, v) oc /(i/max) corrcspouds to the diagonal elements K^i, in the discretized equation ( fTTI) . 
However, for any value of 7, f{ymia^ is extremely small for strong physical reasons: a single 
nucleon has a vanishing probability of carrying the entire momentum of the nucleus. Thus 
the matrix K has very small values along the diagonal and is very close to singular, so this 
solution method fails. 

A standard approach to solving Volterra equations with a kernel vanishing identically 
along the diagonal is to either integrate Eq. (fT2l) by parts, or to differentiate with respect 
to X. The first technique gives an integral equation for the primitive of with kernel 

dK{x,t)/dt, while the second has kernel dK{x,t)/dx and left- hand- side dJ^'^{x)/dx. These 
approaches are still problematic, however, because derivatives of the smearing functions are 
still very small at y = t/max! also, taking derivatives of functions derived from fits to data 
introduces substantial errors. Furthermore, the solution depends on knowing jF"(x) at all 
values of x, while in practice, only data up to x = 1 are available. It is clear that a direct 
solution to Eq. ([9]) is impractical for the particular forms of smearing functions used in this 
model. 



B. Multiplicative solution 



The most widely-used method for extracting spin-averaged structure functions is the 
smearing-factor or multiplicative method 8|]. This is an iterative solution method based on 
the ansatz that the right-hand-side of Eq. iQ can be written as a product of the neutron 
structure function and a "smearing factor" S^{x), 



T'^ix) 



S'^ix) T'^ix^ . 



(13) 
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From a first guess one obtains S^^\x) by smearing J-'"'^^\x) and dividing by 

jF"(°)(x). Dividing jF"(x) by 5""^°-' gives so the result after one iteration is 

One can see from the form of Eq. ( fT4l) that this method is problematic if the smeared struc- 
ture function has zeros in the range of x of interest. The spin- averaged nuclear structure 
functions are positive-definite for < x < 1, so this problem does not arise, and the mul- 
tiplicative method converges quite rapidly for essentially any reasonable choice of jF"(°)(x). 
Even for spin-dependent structure functions, which may have several zeros, the multiplica- 
tive method works fine as long as the zeros of the smeared JF"^'') (x) are very close to the zeros 
of Since the smearing functions are close to 5-functions, this amounts to requiring 

that the zeros of the neutron structure function be very close to the zeros of the nuclear 
structure function. Experimental errors could easily obscure the true location of the zeros 
of the nuclear structure function, though, making a direct application of this method to 
experimental data difficult. 

C. Additive extraction method 

Instead of assuming a multiplicative smearing factor, one can exploit the fact that the 
smearing function / is sharply peaked about y = 1 to formally write 

/(l/,7) = ^^6iy-l) + 5/(1/, 7) , (15) 

where J\f = j^^^^ dy f(^y^ ^) jg the normalization of the smearing function, which for 7 = 1 
is either unity for the unpolarized F2 structure function, or equal to the effective nucleon 
polarization in the nucleus for the spin-dependent gi structure function. The correction 6f 
gives the finite width of the smearing function. The smeared neutron structure function in 
Eq. iQ can then be written 

^-(x) = Af r^ix) + (5/®J^")(x). (16) 

The convolution term in Eq. (fT6!) can thus be treated as a perturbation and the equation 
solved iteratively. Starting from a first guess JF"'^'') (x) one has, after one iteration, 

1 
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^"(x) - (/®^"(°)) (x) . (17) 



Here, there is no danger of divergences due to zeros in the input, as the only division is by 
A/", which is nonzero for all smearing functions f{y,'y)- 

When = xgi, Eq. ([H]) is a system of two equations, whose solution is slightly more 
involved. One notes that the function fu is the most sharply peaked of the smearing 
functions [l^, and hence gives the largest contribution to xgf. Assuming that the fu 
contribution is zero, one can apply Eq. (fTTl) to xgf to obtain gi^^\x), which is substituted 
into the expression for xg2- Subtracting this contribution /21 ® a^fi'" from xg2 and applying 
Eq. (fTTl) to the resulting expression then gives g^^^\x). The new value g2'^^\x) is then 
inserted into the xgf equation and the recursive procedure repeated until convergence is 
achieved. 



D. Analysis of convergence 

As we will show in Sec. IIVI below, the convergence of the additive method is quite fast and 
nearly independent of the initial guess. The reason for this is essentially the sharply peaked 
shape of the smearing function. This can be illustrated by examining the propagation of the 
error on the true function J?^truc(a;) with each iteration i. Starting from a first guess, 

^(0)( 

for the true function, we define !F^^\x) = J-'truc{x) + e^^\x), where e^'^^x) is the difference 
between the first guess and the true result. Tracking this error after i = 1 iteration gives 

g(i) (3,) = g(o) - ^ (/ ® e(°)) (x) . (18) 

Note that if /(y) ^ 6{y — 1), the error for a; < 1 vanishes even after one iteration. In 
fact, since f{y) is sharply peaked at y = 1 (for the unpolarized and diagonal polarized 
distributions), e^^\x) is expected to be quite small for x < 0.8, regardless of J-'^'^\ 

More specifically, the iteration procedure will converge if for successive iterations 
^ |e'^*)(a;)|. Defining a to be the width over which the smearing function f{y) ^ 0, 
from Eq.(j4]) it follows that cr ~ ^Pchai/M <^ 1, where Pchar is the characteristic nucleon mo- 
mentum scale in the deuteron. Then using the generalized mean value theorem for integrals, 
the correction term in Eq. (fTSl) can be written as 

^{f®e){x) = e{x/y.), (19) 

where = 1 + c is a point within the integration interval, with |c| < cr/2. If e(x) is a 
sufficiently smooth function of x, one can expand the right-hand-side of Eq. f|T9l) in a series 

11 



in c, 

e(x/(l + c)) = e(x) - cx e\x) + 0{c^) , (20) 
so that the error after one iteration is e'^^^x) ~ cxe'^^'^' [x] . This then leads to the estimate 







g(0)' 


0" 


e(o)' 


e(0) 


^ cx 


e(o) 


< - 
2 


e(o) 



so that the ratio of errors e'^^^^^'^'' proportional to the width of the smearing function, as 
long as the width is small. Furthermore, because the i = 1 error is given by the derivative 
of e*^''-*, convergence is fastest when the error is smoothest, which will typically be away from 
resonance peaks. 



IV. RESULTS 



In this section we present numerical results which illustrate the features of the extraction 
methods discussed in Sec. IIIII We discuss firstly the unpolarized structure function, 
before considering the more challenging case of the polarized structure function. 



A. Unpolarized structure functions 

Most previous extractions of the structure function have been performed in the deep 
inelastic region, where the structure functions are smooth and monotonic (beyond x ~ 0.3). 
Before applying our extraction procedure to the more challenging resonance region, we first 
test the method on the more familiar case of DIS kinematics. For the input proton and 



neutron structure functions we use the MRST parameterization 251] at = 10 GeV^, 
and simulate the deuteron "data" using the finite-Q^ smearing function fo{y,'~f) from 



Ref. [16 1 . 



The resulting extracted neutron F2 structure function is shown in Fig. [2] using the ad- 
ditive method. Starting from an initial guess of F2^^^ = 0, the extracted curve is almost 
indistinguishable from the input F2 after just a single iteration. The main reason for this 
fast convergence is the fact that the nucleons in the deuteron are weakly bound and have 
small average momentum, which leads to a smearing function fo{y, 7) that is sharply peaked 
around y = 1. Although the precise height and width of the peak may vary slightly for differ- 
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FIG. 2: (Color online) Extraction of the neutron F2 structure function at = 10 GeV^ from 



F2 (dashed) and F2 (dot-dashed) data simulated from the MRST parameterization 



smearing function fo{'j,y) 



251 ] and the 



161 ] using the additive method. The extracted F2 structure function 



after i = 1 iteration (solid) is almost indistinguishable from the input (dotted). 



ent deuteron wave functions, the rapid convergence is a relatively model-independent feature 
of the extraction. 

While the extraction of in the deep inelastic region is straightforward, obtaining 
in the nucleon resonance region, where the cross section is dominated by resonance peaks, is 
more problematic. In fact, to our knowledge such an extraction has not yet been undertaken 
in any quantitative analysis. Even in a system as dilute as the deuteron, the structure of 
nucleon resonances is significantly smeared out by the Fermi motion of the nucleons, so 
that for ~ 1 GeV^ or higher essentially only the A region exhibits any clear resonance 
structure. In heavier nuclei there is very little resonance structure evident at all It is 
not clear a priori therefore to what extent neutron resonance data can be extracted from 
data in which the neutron information is strongly smeared. 

To test the effectiveness of the additive extraction method in the resonance region we use 



as input structure functions from the MAID Unitary Isobar Model 



26|, which is constructed 
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FIG. 3: (Color online) Extraction of the neutron st ructure function at = 1 GeV^ in the 



resonance region using the MAID parameterization 



26| for the input (dotted), and a first guess 



n"{0) 



0. The convergence of the procedure is illustrated by the results for i = 1 (dot-dot-dashed), 



2 (dot-dashed), 5 (dashed) and 30 (solid) iterations. 



to parameterize meson electroproduction data at low W. The convergence of the iteration 
procedure in the resonance region is illustrated in Fig. [3l where we attempt to extract the 
input F2 at Q"^ = 1 GeV^ with an increasing number of iterations. Taking as a first guess 



,n{0) 



0, after i = 1 or 2 iterations the prominent resonant structures are clearly visible, 



although the amplitudes of the resonance peaks is still underestimated. After i = 5 iterations 
the extracted function is very close to the true result, and would in most cases lie within 
experimental uncertainties. Repeating the procedure i = 30 times reproduces the complete 
resonance structures almost exactly. 

The multiplicative method can also be used to extract F2 in the resonance region, as 
illustrated in Fig. |H The starting point for the iteration here is taken to be F^*'^'' = F^, 
and after i = 5 iterations the result is in good agreement with the input function, only 
slightly underestimating the peaks of the resonances. As in the additive method, almost 
perfect agreement can be achieved eventually with further iterations. Note that a direct 
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FIG. 4: (Color online) Comparison of the convergence of the additive and multiplicative methods 
for extraction F2 for i = 2 and 5 iterations. The starting point for the additive method was 



n"{0) 



0, while for the multiplicative method -Fg 



n(0) 



taken from the MAID parameterization 



F2 . The input structure functions were 



2g at Q2 = 1 GeV2. 



comparison of the convergence of the additive and multiplicative methods from Fig. H] is not 
possible since the starting points F^*-'^'* are different. Here we merely illustrate the fact that 
both methods can converge to the true result within a relatively small number of iterations. 
To examine the sensitivity of the extraction to the initial guess F2^^^ , in Fig. [5] we show 



the result after i = 2 iterations for initial guesses = and 



,n(0) 



fit 



F|, using the MAID 



26| at = 1 GeV^ as input. Since the amplitudes of the resonances are significantly 
larger for the proton than for the neutron, the initial guess results in larger amplitudes 
for the extracted neutron F2 for the same number of iterations. On the other hand, because 
the proton and neutron resonance transitions to the A are expected to be equal (since the 
transitions are isovector), the proton initial guess enables the A peak to be reproduced 
extremely well, in contrast to the zero first guess which requires more iterations to produce 
the observed structure. Of course, with sufficiently many iterations the input F2 can be 



reproduced accurately regardless of the initial guess F2 



n(0) 
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FIG. 5: (Color online) Extracted neutron F2 structure function using the additive method with 
i = 2 iterations, starting with initial guesses F^^'^^ = (labeled "n(0) = 0", dashed) and F^^^^ = 
F2 (labeled "n(0) = p", solid). The input neutron (dotted) and proton (dot-dashed) structure 



functions are taken from the MAID parameterization 



261] at Q2 = 1 GeV2. 



In all of the above extractions the full 7- (or Q^-) dependent nucleon smearing function 
/o(?/, 7) has been used when computing the deuteron structure function. While using a 
7-independent smearing function may be a reasonable approximation in the deep inelastic 
region where 7 values are typically close to unity, applying the 7 = 1 smearing function 
to low-Q^, large-x data can lead to errors in the extracted F^, especially in the resonance 



region 28 



The importance of using the correct smearing function is illustrated in Fig. [HI where 
we show the extracted neutron structure function after i = 10 iterations. The result 
using the full, 7-dependent smearing function is very close to the input. On the other hand, 
with the Q^-independent, 7 = 1 smearing function the iteration does not converge to the 
correct solution. In particular, while a resonance bump is visible in the A region, it has 
the incorrect strength; the second resonance region displays a trough where there should be 
a peak; and the third resonance region appears to have no structure at all. Increasing the 
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FIG. 6: (Color online) Extracted neutron structure function using the additive method after 
i = 10 iterations with the full 7-dependent smearing function (solid) and with the 7=1 approx- 
imation (dashed), compared with the input neutron (dotted) structure function from the MAID 



parameterization 



26!] at Q2 = 1 GeV2. 



number of iterations for the 7-dependent smearing function leads to ever closer convergence 
to the input F^. For the 7 = 1 smearing function, the result does not change qualitatively 
with further iterations, however significant noise develops over much of the x range. 

These features arise from the mismatch between the smearing functions used to compute 
the deuteron and those used to perform the extraction. Of course, had the deuteron 
structure function been simulated with the 7 = 1 smearing function, the extraction with 
the same function would return the same input F2 as in Fig. [6l However, this comparison 
demonstrates the sensitivity of the extraction to the dependence of the smearing function, 
and highlights the importance of using a smearing function with the correct dependence 



when analyzing actual data 
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FIG. 7: (Color online) Deuteron xgf and structure functions simulated from the MAID pa- 
rameterization of the proton and neutron gi^2 at = 1 GeV^ and the smearing functions fij, 



i,j = 1,2, from Ref. 



181 ]. The full results for xgf (solid) and (dashed) are compared with the 



diagonal approximations using fn only (dotted) and /22 only (dot-dashed), respectively. 
B. Polarized structure functions 

At finite the complete expressions for the gf and g2 nuclear structure functions in 
Eq. ([3]) represent a coupled set of equations involving contributions from both the and 
g2 structure functions of the nucleon [l^. While the diagonal /n and f22 smearing functions 
dominate for most kinematics, the off-diagonal fu and /21 contributions could be important 
at low values of Q^. Furthermore, the g^ contribution to g2 survives even in the Bjorken 
limit. 

As described in Sec. IIII CI above, one can solve such a system of equations by simultane- 
ously iterating both g^ and gl^, given known (or simulated) proton and deuteron data. Such 
a procedure will necessarily be slower and require more iterations, but is stable and will in 
principle converge to the correct solutions. 

In practice, however, for the kinematics discussed here, namely ~ 1-10 GeV^, the off- 
diagonal contributions are rather small. This can be seen in Fig. [7] where we show the xgf 
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and structure functions simulated from the MAID gl'2 parameterizations 26| at = 
1 GeV^, using the smearing functions fij{y, 7), i,j = 1, 2, from Ref. 13]. The results with the 
diagonal terms only (/n for gf and for (72); ^-re very close to the full results which include 
both diagonal and off-diagonal contributions. With the precision achievable in current and 
near-term future experiments, the diagonal approximation to the gf2 structure functions 
should therefore provide a reliable framework in which to extract neutron structure functions, 
and in the following analysis we consider only the diagonal contributions. Furthermore, 
since the shape of g2 is qualitatively similar to that of gi (generally g2 has the opposite sign 
compared with gi), we shall focus on the gi structure function as representative of the effects 
of extracting spin-dependent neutron structure functions in the resonance region. 

As we saw in the previous section, both the additive and multiplicative methods yield 
reliable results for extracted neutron structure functions, in both the deep inelastic and 
resonance regions, as long as the structure functions are free of zeros. For polarized scattering 
the gi and g2 structure functions are no longer positive-definite, so that taking ratios of 
smeared to unsmeared functions can in principle lead to singularities during the extraction. 

This does not necessarily render the multiplicative method completely impractical for 
extracting polarized structure functions, however. Numerically, for a given iteration where 
the structure function is close to (but not exactly at) its zero, the smearing factor S"" will be 
very large. For the next step in the iteration this large contribution will be damped by the 
corresponding small value of the structure function, making the result finite. On the other 
hand, precisely how (and whether) this cancellation occurs in practice will be determined 
by the shapes of the input structure functions and smearing functions, and a priori it is not 
clear whether an extracted non-positive definite structure function will be well-behaved for 
a particular extraction. 

To illustrate the extraction of spin-dependent structure functions we first consider the gi 
structure function in the DIS region in Fig. [SI The input proton xgi and neutron xg^ data 
are taken from the leading twist parameterization in Ref. 30| at = 10 GeV^, with the 
deuteron xgf simulated using the smearing function fii{y, 7) from Ref. 18|. With a starting 
point of xgi^^^ = 0, the extracted neutron structure function after a single iteration using 
the additive method is essentially indistinguishable from the input. As for the unpolarized 
F2 structure function in the DIS region in Fig. [2], this feature reflects the narrow width of 
the smearing function /ii(?/,7) around y = 1. 
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FIG. 8: (Color online) Extraction of the neutron xg\ structure function from proton xg\ (dot- 
dashed) and deuteron xg\ (dot-dot-dashed) "data" simulated from the leading twist parameter- 



ization 



30( 1 at = 10 GeV^ and the smearing function /ii(7,y) 



161]. The input (dotted) and 



extracted xg^ functions, for both the additive (solid) and multiplicative (dashed) methods, are 
almost indistinguishable. 

For the multiplicative method the initial guess is taken to be xg^^'^ = xg^, and after 
one iteration the extracted neutron structure function is also very close to the input. In 
particular, even though the ratio of smeared to unsmeared g^ structure functions is singular 
at a; ~ 0.5, the extracted function is nevertheless continuous in this region. The marginally 
slower convergence here compared with the additive case reflects the different starting inputs 
for g^, which for the multiplicative method is further from the true result than for the 
additive. 

While both the additive and multiplicative methods appear to be effective in extracting 
the spin-dependent neutron structure function in the DIS region, their utility in the nucleon 
resonance region, where the xgi exhibits considerably more structure, is compared in Fig. [HI 
using the MAID parameterization 26| at = 1 GeV^. The most striking feature of the 
extracted neutron xg'i is the discontinuities near the zeros of the input function for the 
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FIG. 9: (Color online) Comparison of the convergence of the additive and multiplicative methods 
for extraction xg" for i = 2 and 5 iterations. The starting point for the additive method was 
xgl^'^^ = 0, while for the multiplicative method xg^^'^^ = xg\. The input structure functions were 
taken from the MAID fit [201 at = 1 GeV^. 

multiplicative method, which arise from the singularities in the smearing factor S"". On the 
other hand, no such singularities appear for the additive method and the extracted functions 
are smooth and continuous over the entire range of x. 

The convergence of the extraction for the additive method is illustrated in Fig. [10], where 
after only five iterations the extracted xQi displays all of the prominent features of the 
A peak and the higher resonance regions. After z = 30 iterations the input function is 
reproduced almost exactly. 

The sensitivity of the extraction to the initial guess for xg^ is illustrated in Fig. [HI where 
the results after i = 2 iterations are compared for starting values xg^^^^ = and xgi^^^ = xg\. 
As in the case of the unpolarized F2 structure function, the xg\ initial guess gives amplitudes 
that are larger than for the zero initial guess after the same number of iterations. Since the 
input proton and neutron structure functions are similar in the A region, the iteration of xgi 
converges on the A peak more rapidly for the xg\ starting point than for the zero first guess. 
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FIG. 10: (Color online) Convergence of the extracted neutron xg" structure function for z = 1 
(dot-dot-dashed), 2 (dot-dashed), 5 (dashed) and 30 (solid) iterations, using the MAID resonance 



fit 



2^ at = 1 GeV^ as input (dotted), with a first guess xg^^^'' = 



On the other hand, because the second resonance peak for the proton is significantly larger 
than for the neutron, convergence on this is faster for the xgi^'^^ = starting value. Again, 
with sufficiently many iterations the input xg^ can be accurately reproduced independently 
of the starting point. 

The importance of using the correct dependence in the smearing function fu{y,l) is 
highlighted in Fig. [121 where the extracted xQi neutron structure function is shown after 
i = 10 iterations. While the full, 7-dependent smearing function yields an almost exact 
reconstruction of the input structure function, the result using the 7 = 1 smearing function 
bears little resemblance to the true xQi- Most notably, the height of the A peak is signifi- 
cantly underestimated, and the position of the second resonance peak does not correspond 
to the correct value. As for the F2 structure function in Fig. El these features arise from the 
mismatch between the smearing functions used to compute the deuteron xgf and those used 
to perform the extraction of xgi. They clearly demonstrate that it is vital to use the cor- 
rect dependence in the smearing function when analyzing data in the nucleon resonance 
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FIG. 11: (Color online) Extracted neutron xg^ structure function using the additive method 



with i = 2 iterations, starting with initial guesses xg^ 



n(0) 



(labeled "n(0) = 0", dashed) and 



n(0) 



xg\ (labeled "n(0) = p", solid). The input neutron (dotted) and proton (dot-dashed) 



structure functions are taken from the MAID fit 



26| at Q2 = 1 GeV^. 



3l| 



region, especially at low Q"^ and large x 

Analysis of actual gf (and F^) data to extract the free neutron structure functions will be 
discussed in a forthcoming publication 29|]. However, one can anticipate how the neutron 
xg^ structure function can be extracted from actual proton and deuteron data, together 
with error bars, by a simple illustration. 

In Fig. [13] we show the proton and deuteron structure functions simulated from the 
leading twist parameterization 30|] at = 10 GeV^, with the error bars derived from 



the uncertainties on the xgi and xgi structure functions given in Ref. 30|]. The deuteron 
structure function was simulated by varying each point of the proton and neutron input by 
a Gaussian of width given by the error bar, which were then smeared with the momentum 
distribution fu and added to get a "trial" xgf. This procedure was repeated for 50 trials, 
after which the average and standard deviation of each point was taken to obtain the xgf 
curves and error bars. 
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FIG. 12: (Color online) Extracted neutron xg^ structure function using the additive method 
after i = 10 iterations with the full 7-dependent smearing function (solid) and with the 7=1 
approximation (dashed), compared with the input neutron (dotted) structure function from the 



MAID parameterization 



261] at Q2 = 1 GeV2. 



To extract xr?", one can assume that the only errors that contribute are those from the 
deuteron. In practice, errors on xgf are much larger than those on xg^, and smearing xg\ 
renders the proton errors negligible compared to the deuteron errors. As before, each point 
of xgf was varied by the error bars, from which the smeared xg^ was then subtracted with 
no errors, and the extraction performed to obtain a "trial" xg^. This was repeated for 50 
trials and the average and standard deviation computed as before. 

The extracted neutron data points in Fig. [T3] are found to be in excellent agreement 
with the input xg^ structure function. The errors on the extracted neutron function after 
one iteration are of the same order of magnitude as those on the deuteron. Note that 
the extracted error bars are considerably smaller than the original error bars, indicated by 
the shaded band around the input xg"^, which is mostly due to the fact that, as a sum of 
smeared functions, the simulated xgf has artificially small errors. For real data, errors on 
the deuteron and proton structure functions are given, and neglecting the errors on xg\ with 
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FIG. 13: (Color online) Extraction of the neutron xg^ structure function (circles) from proton xg^ 



sc 



uares) and deuteron xgf (diamonds) "data" simulated from the leading twist parameterization 



30| at = 10 GeV^ and the smearing function /ii(7,y) 
the uncertainties on the structure functions given in Ref. 



16l |. The error bars are derived from 



30]. The input xg'l structure function 



(solid) is given as reference, with uncertainties indicated by the shaded band, 
respect to xgf is a very reasonable assumption. 

V. CONCLUSIONS 

In this paper, we have presented a new method which allows the reliable extraction of 
neutron structure functions, both spin-averaged and spin-dependent, over a wide range of 
Q^. We have compared the new (additive) method to the existing (multiplicative) extraction 
method, and found that the performance of both methods is very similar for the extraction 
of F2, while the additive method is free of the singularities that develop when attempting 
to extract xg^ using the multiplicative method. Moreover, the speed of convergence of the 
additive method is nearly independent of the initial guess, and in most reliable 
extraction is achieved after i = 5 iterations. Finally, the extraction of xg"^ including error 
bars, shown in Fig. [131 illustrates both that errors on the extracted function can be reliably 
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estimated, and that the performance of the method is not overly sensitive to perturbations 
in the input. 

The tests of the additive method on models of resonance-region structure functions show 
that the general shape of the curve is reproduced after only one iteration, but that further 
iterations are necessary to accurately extract the magnitude of the resonance peaks. The 
dependence of the method on the initial guess is evident in the sense that the regions where 
convergence is slowest are the regions where the initial guess is farthest from the actual 
magnitude of the resonance peaks. Thus, one can reduce the number of iterations needed 
with an educated guess about the shape of the neutron function. For example, since the 
isovector transition to the A gives identical proton and neutron structure functions for the 
resonant part of the A, a good first guess for the neutron would always be the proton 
structure function in the A region. To ensure that the extracted neutron structure function 
is in fact correct in the context of the smearing-function model, one can smear the extracted 
structure function and add to the smeared proton structure function to compare with the 
deuteron data. 

Despite extensive experiments on light nuclear targets, the neutron remains something 
of a mystery. The same observables which can be directly measured for the proton must be 
inferred for the neutron, because its instability outside of the nucleus makes neutron targets 
impossible. Previously, the low statistics and large errors from experiments designed to 
measure neutron observables limited the accuracy of measured neutron structure functions 
far more than using a simplified model of the nucleus to perform the extraction. The 
situation has changed with recent experiments at JLab, and now accounting for nuclear 
corrections in neutron structure function extraction procedures is essential to obtain an 
accurate representation of the neutron structure functions, especially in the resonance region. 
In particular, we have shown that ignoring finite-Q^ corrections to nuclear structure functions 
leads to an extracted neutron structure function which may bear little resemblance to the 
true shape. To assess quark-hadron duality for the neutron to the same extent that it has 
been verified for the proton, detailed knowledge of all neutron structure functions in all 
kinematic regimes is needed. The method presented in this paper, when applied to the most 
recent JLab data, will be a first step in that direction. 
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